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Abstract 



We suggest a simple physical picture for the diffractive parton distributions 
that appear in diffractive deeply inelastic scattering. In this picture, partons 
impinging on the proton can have any transverse separation, but only when 
the separation is small can they penetrate the proton without breaking it up. 
By comparing the predictions from this picture with the diffractive data from 
HERA, we determine rough values for the small separations that dominate 
the diffraction process. 
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Typeset using REVT^ 



Diffractive deeply inelastic scattering (DIS) has been extensively studied at HERA [|l|. 
For the theoretical description of this process, one can usefully factor the perturbative short 
distance physics from the non-perturbative long distance physics. In such a description, the 
diffractive part is long distance and is described by diffractive parton distributions. The 
measured diffractive structure function -F^'^ is a convolution of these distributions with the 
usual hard scattering coefficients known from inclusive DIS. This is described in detail in 
Ref. and the references cited there. 

One of the results of Ref. is that the diffractive gluon distribution, although it con- 
tributes to diffractive DIS only at the next-to- leading order in a^, is extremely important in 
diffractive DIS because it is very large. The magnitude of its ratio to the quark distribution 
is found to be governed by the color factor C\{N^ — l)/{CpNc) = 27/2. A very large gluon 
distribution makes diffraction an important process to probe the onset of parton satura- 
tion 1^-^] and the limits of applicability of perturbative evolution 0. Recent experimental 
observations of diffractive jet production [|^, which are especially sensitive measurements of 
the gluon distribution, confirm the existence of a very large diffractive gluon component. In 
this paper, we build on the results of Ref. to develop a physical picture for the diffractive 
gluon distribution based on the analysis of absorption (color opacity) and scattering (color 
transparency) effects in the diffraction process. 

The diffractive gluon distribution is a product of matrix elements of the gluon measure- 
ment operator 



Ga{0,y~,Oy^= Vexp(-tg dx' A^{0,y- ,0%) G,{0,y~,O)+^ 



(1) 



ab 



where, in standard notation, G is the gluon field strength, A the vector potential, t the 
generators of the adjoint representation of SU{3), and V denotes path ordering of the ex- 
ponential. As explained in Ref. p|, when viewed in the proton rest frame this operator 
creates a state consisting of partons that have large momentum in the minus direction 
(/c^ = (fc° ± k^)/\/2). These partons travel a long distance, then impinge on the proton 
and scatter. The operator in Eq. (|l]) contains an eikonal line operator that is capable of 
absorbing gluons from the color field of the proton and is equivalent to a color octet parton 
moving with infinite minus-momentum. Note the correspondence with the space-time de- 
scription of diffractive DIS [p|-p!l[], in which a current operator J'^(0, y~, 0) creates a state of 
large momentum partons that impinge on the proton and scatter. In fact, in the diffractive 
gluon distribution, the color octet eikonal operator is a stand-in for a quark-antiquark pair 
that was produced by the current but is too compact to be resolved by the color field of the 
proton. 

The partonic system created by the gluon operator can have any transverse size when 
it reaches the proton. We propose here that, if the partons reach the proton at a large 
separation, they cannot get through the proton without breaking it up. Thus diffractive 
deeply inelastic scattering results mainly from the incoming partons being close together and 
interacting with the color field of the proton only through their color dipole moment. This 
would seem to make diffraction quite unlikely, but the calculation shows that the partonic 
wave function is large at small transverse separations between the partons. Because of the 
short distance behavior of the perturbative wave function, diffraction turns out to be not at 
all unlikely. 
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In the following, we first restate the result p[ for the diffractive gluon distribution. 
Then we formulate the picture for the interaction of the partonic system created by the 
gluon operator with the proton's color field. We use this in conjunction with the short 
distance factorization of the hard scattering to calculate the diffractive structure function 
F2^^. Finally we compare the results to experiment and extract a rough value for how close 
together the partons have to be in order to make it through the proton without breaking it 
apart. 

In the notation of Fig. 1, the expression for the diffractive gluon distribution is 

dx^dt 167r3 ^{l-P)xl2^^^,J (2vr)2 ^ ^ 

with 

= J J^$^ J dy-e'^^-^iy-{pA',SAr,k,s,a\Ga{0,y',oy^pA,SA) . (3) 

The expression in Eqs. (0) and represents a perturbative approximation in which the 
high momentum gluon created by the operator G travels to the final state without splitting 
into any further high momentum partons, so that there is a single final state gluon with 
momentum k in addition to the final state proton. If one worked to higher orders of pertur- 
bation theory, then there could be partonic splitting and joining and the final state could 
contain multiple high momentum partons. We integrate over the transverse momentum k 
of the final state gluon; its plus momentum is set to k'^ = {1 — f3)xpp\ so that the total 
plus-momentum removed from the proton is Xpp\. Its color h must be 6 = a since the initial 
and diffracted proton states are color singlets. 




FIG. 1. Structure of the diffractive gluon distribution. The top part is treated perturbatively. 
Any number of soft gluons are exchanged with the hadron matrix element represented on the 
bottom. 

At large 1/xp an eikonal approximation applies. There is a wave function ■?/', describing 
the state of the fast moving gluon as it reaches the proton. This gluon passes through the 
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proton at transverse position A while the eikonal operator in Eq. ([l|) is located at transverse 
position 0. The result, taken after a little manipulation from Eq. (3.38) of Ref. f^, is 



iV2-l J (27r)2p+ y A2 i (27r)2'' di^ di^^ {P,k,k + z) 
x{pa',sa'\Tt\1-F{A)F^{0)]\Pa,sa) • (4) 



Here the operators F are color octet eikonal line operators: 

F(z) = V exp l^-ig £j dz~A^{0, , z) t,| . (5) 
The wave function for the gluon state is 

We choose here to differentiate ip twice with respect to £. This has the effect of eliminating 
apparent ultraviolet divergences from the integration over k in the formulas that follow. The 
derivatives with respect to i are compensated by a factor (— l/A^) multiplying the matrix 
element of the eikonal operators. The factor 1 / A^ does not itself produce a short distance 
divergence because both the matrix element and its first derivative with respect to A vanish 
at A = 0. See Sec. 4.1 of Ref. for a discussion of the ultraviolet behavior in the minimal 
case of two gluon exchange; exchanging more gluons does not make the ultraviolet behavior 
worse. 

Consider now the matrix element of the eikonal operators between the initial and final 
proton states. The initial state carries zero transverse momentum. Let — q be the transverse 
momentum carried by the final state (q^ ~ —t) . We may first expand this state in terms of 
eigenstates of transverse position with transverse position — b as 

W, sa'I = {Pa', -q, ^A'\= J d% e''^ "^ ^{p^,, -b, sa'\ , (7) 

where the subscript x on a state vector denotes that it is a position eigenstate. We then 
apply a translation to put the proton at transverse position and let the eikonal operators 
be at transverse positions b and b + A. Taking into account that the eikonal operators do 
not transfer any plus momentum, we define a non-perturbative function S as 

(27r) 2pt 5ip+ - pt) H(b + A, b; sa', sa) 

^ 0, SA'lTr [l - F(b + A)F^(b)l 0, sa) . (8) 



iV2-l 



We use Eqs. (0) and (D in Eqs. (|) and (|) to give 

dx^dt ~J A2 ^,2dbdbe l[x^,(J,/X,/X) 
x^Y. ^*ih' + A', h' ISA', SA)E{h + A,h;SA', sa) (9) 
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where 



J(a:,,/3,A,A') 



1 



16 7r3 p{l-p)x^ 



d d 

— r — zi^^%p,'k,k + e') 



2 / (27r)2 

d_d_ 



dH ^_,t^ 



(27r)2 



dH' 



(10) 



Eq. (0) provides a representation for the diffractive gluon distribution in a physically 
appealing form. The (3 dependence is entirely contained in the factor X, which is constructed 
from the perturbative wave functions if). There are integrals over the global impact parameter 
of the collision (b for the amplitude and b' for the complex conjugate). The factor I does 
not depend on b or b'; these parameters merely tell where the proton is relative to the 
incoming partons. There are also integrals over the parton transverse separation (A for the 
amplitude and A' for its complex conjugate). The perturbative function X/(A^ A'^) gives 
the amplitudes for the partons to have separations A and A'; then the functions S describe 
the effect of the proton's soft color field. 

If we integrate the diffractive gluon distribution (P) over Xp at fixed x (or equivalently over 
/?), this integration produces a logarithmic divergence coming from the region (3 = x/xp <^1. 
In the case of the quark distribution there is no such divergence 0. Then by integrating 
over (3 one obtains a formula with only one separation parameter A. Formulas of this kind 
are given, e.g., in Refs. P,P,p!l|] . For the discussion that follows we are interested in the full 
13 spectrum, for which formulas with two separation parameters, A and A', apply. 



If we integrate the diffractive gluon distribution (^ over t = — we obtain 



df;;fMxp,fi) 



An 



d'^A d^A' 



^d^b J(xr,/5,A,A') 



1 



X- ^ H*(b + A',b;sA',SA)S(b + A,b;s^/,s^). 

SA,SaI 



(11) 



Now there is only one impact parameter, b, but there are still two parton separation vari- 
ables, A and A'. 



We now turn to the hadronic matrix element S. If, as in Ref. [jl2l, the proton were 
replaced by a small meson made of heavy quarks, then perturbation theory would be appli- 
cable to compute S. But for a large state, such as a proton, we need to rely on a combination 
of perturbative and nonperturbative ideas to construct an approximation for S. 

We first recognize that S is an absorption amplitude. Imagine that two high momentum 
gluons, represented by eikonal operators, impinge upon the proton at transverse positions 
b and b + A. Suppose that the gluons completely miss the proton. Then each of the two 
operators F in 1 - F(A)F'^(0) can be replaced by unit operators so that S = 0. Now, 
suppose that one or both of b and b + A is inside the proton, with no particular relation 
between these positions. Then we expect that the amplitude for the proton to remain a 
ground state proton after the passage of the gluon(s) is approximately zero, 
so that E ^ 1 X 5 



Thus a simple model for S might be 



^(b + A,b) S, 



SASa' 



(12) 
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where ^ is a smooth function such that ^~Oif|b + A| ^ Rp and |b| ^ and ^ 1 if 
either |b + A| <^ i?p or |b| Rp, where Rp is the radius of the proton. 

An important exception arises if A is small. Then F{h)^ is nearly the inverse of F(b+A), 
so that S f» 0. To take this color transparency effect into account, we modify our model to 

^(b + A, b) 5,,,^, [1 - exp(-/s:2AV2)]. (13) 

Then S ^ for small A but for large A we retain the absorption physics explained above. 
Here l/ng is a parameter that measures how close together the gluons have to be before 
they can traverse the proton without breaking it up. That is, l/ng is the largest size for a 
color octet dipole for which the proton is transparent. We hypothesize that l/ng is small: 

1/ Kg <C Rp. 

In Eq. (|13D we have chosen a gaussian function for reasons of simplicity. Any function 
that vanishes for k^A ^ 1 and tends to 1 for k^A ^ 1 could be used in place of [1 — 
exp(-«:2AV2)]. 

We now note that the function X is singular when A and A' become small. In fact, 
straightforward dimensional analysis shows that for |A| ~ |A'|, X ~ const. /(|A| |A'|). 
For this reason, the important contributions to the integral over A comes from small A, 
|A| ~ I/kj <^ Rp. In this circumstance, we can set A to zero in ^(b + A, b), giving 

H(b + A, b; SA', sa) = A{h, b) 6s^s^, [1 - exp{-KlAy2)]. (14) 

This is our model. The proton can survive, S ^ 0, if either b is large or A is small. The 
model incorporates saturation P,^,|Tl|: as |A| increases, the absorption amplitude increases, 
but it never gets to be larger than 1. Furthermore, the important transverse separations are 
small 

We now use Eq. (|Hp in Eq. ([TI|). The integral in b over the absorption function A must 
give a quantity of order R"^: 

J d% \A{h,h)\' = Rl , (15) 

where the precise coefficient A_4 depends on the detailed form of the function A. In the 
numerical evaluations that follow we take -Rp to be the charge radius of the proton, Rp = 
0.88 fm |[T^, and we take = 1. Using Eq. (|15|) and Eq. (^) for X in Eq. ([Tl|), we see that 



dfSfA _ {N^-1))^aRI^I 



9{P) , (16) 



where g{(3) is a function of (3 that we determine by calculating the integrals numerically in 
a suitable momentum space representation. The behavior in Kg is perhaps surprising - the 
larger k^, the larger the cross section. But it reflects an essential feature of the diffraction 
process: because of the singular nature of the wave functions the amplitude is dominated 
by distances of order I/k^. 

Notice that the result (p!6| ) has an Xp dependence This results from our simple 

treatment of the soft exchanged gluons, which we treat as having rapidities near zero in the 
proton rest frame. In a more sophisticated treatment, one would have a BFKL ladder of 
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exchanged gluons with rapidities becoming large and negative along the ladder. This would 
lead to a dependence with 6 proportional to a^. We do not pursue such a treatment 

here. However, in order to compare with experimental data, we suppose that the simple 
treatment is approximately correct for Xp of order xq, where Xq is small but not too small, 
and that there are corrections for Xp <^ xq- Thus we multiply our result by Xq^ /x"^^. We use 



the experimental value 5 = 0.15 and take xq = 0.01. This gives 

fdiS 

Ig/A _ v- c ^^^^ _ ^^j^ 



A treatment analogous to that described for the gluon may now be applied to the quark 
distribution, using the appropriate quark operators as described in 0. Similarly to Eq. (p!6[), 
the answer has the form 



dxp X] 



■F 



The function is different than in the gluon case due to the different lightcone wave 
functions ip. The color factor is different. Finally, the scale Kg is different because color 
triplet eikonal operators now replace color octet operators in the functions S. From the 
color structure of these operators, we expect Kg and Kg to be roughly in the relation Kg/ Kg ^ 
Ca/Cf- Note that, together with the color factors explicitly shown in Eqs. ([T7| ) and (|18|) 



this gives back the factor Cl{N^ - l)/{ClNc) = 27/2 of Ref. || in the ratio between the 
gluon and the quark distributions. Note also that, since Kg < Kg, the typical separation 
A in the case of quarks will be larger than in the case of gluons, so that the perturbative 
evaluation that we have used may be less reliable in the quark case. 

The treatment discussed so far does not contain scaling violation. Having worked in a 
short- distance factorization framework, we may include the graphs that give scaling viola- 
tion systematically by using next-to-leading-order hard scattering coefficients and DGLAP 
evolution equations, as in Ref. We then compute predictions for the diffractive structure 
function F^^. 

What, then, should be the starting scale /io for evolution? Evidently, /zq should be 
roughly comparable to the physical scales Kg and Kg. To go further, consider the term 
Pgg ® fg in the evolution equation for the diffractive quark distribution at scale /xq- The 
quark and antiquark have a transverse separation of order l/fiQ. The quark-antiquark system 
is separated from the gluon by a larger distance, which, for us, is of order I /Kg. Thus l//io 
should be somewhat smaller than 1/Kg. 

Our results for F^'*^ depend on Kg, Kg, and fiQ. Approximate values for these scales 
may be extracted by comparing the predictions with experiment. Note that a very precise 
determination is difficult because of the uncertainty in the normalization of the model, 
which arises from the uncertainty in the factors and Xq^ in Eqs. (|17D and ([l8|). This 



normalization uncertainty affects the determination of Kg directly. The parameter /xq and 
the ratio Kg/ Kg influence the shape of -^2^'''^ in (3 and Q^, so their determination is less directly 
affected. 

In Fig. 2 we compare the results for -F^'^ with the data of Ref. [|15[ by plotting data 



minus theory over theory, where the theory in the two plots is the one corresponding to 
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FIG. 2. The data for ^2^'^ from Ref. compared with the theory results described in the 
text for two different sets of scale parameters. 



the two displayed sets of values for l//io, l/ft^c,, l/ft^g- Recall that, for the self-consistency 
of the physical picture that we have discussed, we expect these distance scales to be small 
compared to Rp. For the case of Fig. 2a, we find that the agreement between theory and 
data is reasonable, given the simplicity of the model. The size l//io of the quark-antiquark 
pair here is smallest and the ratio Kg/ of the sizes of the color-triplet to the color-octet 
systems comes out to be in rough agreement with the expectation Ca/Cf mentioned above. 
In Fig. 2b we increase the value of l//io and then adjust l/i^g and l/nq so as to describe the 
data as well as possible. We observe that the data seem to disfavor these larger distance 
scales: the shape of the data, in particular, is not well reproduced in this case. 

It is now worth looking back at Fig. 1 and the structure of Eq. (|Tl|). The function 
S, representing the lower part of Fig. 1 and including the exchange of any number of soft 
gluons, is nonperturbative. We have posited a description for S in terms of an absorption 
function and a color transparency factor. The upper part of Fig. 1 is also, in general, 
nonperturbative, because it still contains arbitrarily large parton separations |A|. Indeed, 
in the upper part of Fig. 1 one would expect to have not just two but several fast partons 
impinging on the proton, each with a large transverse separation with respect to the others. 
But if the separations are large, then absorption dominates the hadronic matrix element, so 
that central collisions are unlikely to leave the proton intact. Thus, the parton transverse 
separations must be smaller than the color transparency scale. If this scale is small enough, a 
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perturbation expansion for the upper part of Fig. 1 is justified. At lowest perturbative order 
there are only two partons (including the eikonal line) and we obtain the calculated function 
X. At higher perturbative orders there can be multiple partons, all close together; these are 
the configurations that we include through radiative corrections to the hard scattering and 
evolution. 

The results of Fig. 2 provide a justification for this expansion and lend support to the 
physical picture of diffractive DIS as being due to color-dipole parton systems penetrating the 
proton at small transverse separations. In particular, a partonic system in a configuration 
of a color octet dipole, which corresponds to the diffractive gluon distribution, must be 
especially small in transverse size. This contributes to the diffractive gluon distribution 
being especially large, as is consistent with the data. The picture described in this paper 
provides a concrete model for the mechanism through which the dominance of small size 



states P,p!2| may be realized in diffractive DIS at HERA. More and better data will be 
valuable to investigate this further and decide whether or not the indications of Fig. 2 are 
confirmed. 
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